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Introduction 


Bently Rotor Dynamics Research Corporation (BRDRC) 
has been conducting research on uid flow-generated rotor 
instabilities in fluid handling rotating machines for over 
seven years. The results are reported in a series of publica- 
tions (see References). 

In this article, a summary of the results is presented in a 
simplified form. Noted here are the bottom line, qualitative 
descriptions of physical phenomena, and approximate 
analytical expressions describing these phenomena. 

What is rotor instability? 

The situation when the shalt of a rotating machine exhibits 
high lateral vibrations which are not correlated to unbalance 
or other external periodic forces or constant radial forces 
(such as gravity) is informally referred to as “rotor 
instability’. These vibrations usually have subsynchronous 
frequencies. 

A more formal definition is that the rotor exhibits seli-ex- 
cited vibrations induced by some internal mechanism 
which, in most cases, transfers rotational energy into shaft 
lateral vibrations. In particular, Nuid flow fincluding gas and 
steam) can play the role of such an energy convertor. 
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Fluid circumferential flow — the main 
contributor to rotor instability 

The physical phenomena can be briefly described as follows 
A shalt rotating in an enclosed fuid environment (such as in 
bearings, seals, or a stator/case) drags the fluid into rotative 
motion, Usually the fluid flow is three-dimensional (spa- 
tial). However, the circumferential component (most often 
shaft rotation generated) may appear quite significant, inde- 
pendently from the other components of fluid motion (racial 
and axial). The circumferential flow generates the dynamic 
effect of fluid dynamic rotating force which, in turn (or 
rather in a feedback loop), drags the rotor into lateral 
vibrations. 


Main characteristics of circumferential flow in 
bearings and seals 


For the purpose of analyzing and predicting rotor instability, 
the circumferential flow in bearings and seals is represented 
by two main factors: the fluid circumferential average 
velocity ralio, A (lambda), and Aud film radial stiffness, 
Ay. Both of these factors are functions of many parameters. 
The most important among these parameters is, however, the 
shaft eccentricity, ie. shaft displacement from concentric 
position closer to the wall of a bearing or seal, The radial 
stiffness Kpg increases with shaft eccentricity (Fig. 1). The 
fuid average velocity ratio decreases (slightly) when the 
shaft operates at low eccentricity and decreases dramatically 
when eccentricity becomes higher (Pig. 2). At high cecen- 
tricity, the shaft approaches the stationary wall and “cuts 
aft” the circumferential flow, While this flow situation 
mainly occurs in radial clearances, as in seals and bearings 
it may also occur in an axial clearance area (such as associ- 
ated with balance pistons). 
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Fluid film radial stiffness as a function of shaft eccentricity: c 
is bearing radial clearance. Similar relationship holds true 


for seals. 
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FLUID AVERAGE VELOCITY RATIO, A 


Typical fluid circumferential average velocity ratio as a func- 


tion of shaft eccentricity inside a bearing or seal. 


Fluid film radial stiffness A p as a nonlinear function of shaft 
cceentricity represents a well recognized characteristic of 
bearings and seals. The fluid circumferential average veloc- 
ity ratio Ais relatively new. It has been identified at BEDRE 
throughout extensive experimental testing [Ref. l-16]. 

The more formal definition is as follows: Fluid circumferen- 
tial average velocity AQ induced by shaft rotation with 
rotative speed Q is the angular velocity at which the fuid 
Now-generated force rotates. 

For any bearing or seal, both functions Xp and A can be 
obtained experimentally [Ref. 4-16] and/or analytically from 
Nuid dynamic theory equations [Ref, 23]. 
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5. Whirl and whip 

Fluid motion-induced self-excited lateral vibrations of shafts 

are known as “oil whirl/whip.’ “steam whip.’ “acrocdy- 
namic whip.” or simply “rotor instability’. All these vibra- 
tons actually belong in the same category. They are 
characterized by forward orbital motion fin the direction of 
shall rotation) and subsynchronous frequency, which is eli- 
ther nearly constant (independent of rotative speed) or is a 
fraction of the rotative frequency, maintaining this fraction 
constant as the rotor speed varies, and when other parame- 
ters of the system arc constant. This fraction is very close in 
value to the fluid circumferential average velocity ratio dis- 
cussed above. The first situation is characteristic for “whip,” 
the second for “whirl”. The constant frequency of whip is 
usually close to the rotor natural frequency, most often cor- 
responding to the first lateral mode (it could also be of a 
higher rotor mode) | Ref. 25]. 
There is a very smooth transition fram whirl to whip wath 
increasing rotative speed (Fig, 3a). Sometimes, when A is 
close to 0.5 and the threshold of stability is higher than dou- 
ble the first natural frequency of the rotor first lateral 
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Figure 3 
(a) Spectrum cascade of the vibrational response during 
startup of a rotor supported in one rigid and one oil-lubricated 
bearing. (b) Spectrum cascade of the vibrational response 
during startup of a rigidly supported rotor with a mid-span 
seal, 
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mode, the whip occurs without being preceded by whirl 
(Fig, 3b), At higher rotative speed, the whip may disappear 
and the rotor stabilizes. Sometimes whip “jumps” into a 
whirl of the second mode (‘jump relates to amplitude and 
frequency) (Fie. 4). Sometimes whip and whirl of the seg- 
ond mide exist simultaneously (Fig. 5), 


6. Threshold of stability 


Threshold of stability is the rotor speed at which self-excited 
vibrations start (Pigs. 3, 4, and 5). Very informally speak- 
ing, at the threshold of stability, damping-lype forces disap- 
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Spectrum cascade of the vibrational response during startup 
of a rotor supported in an oil-lubricated bearing indicating 
the “jump” of seli-excited vibrations from whip to the second 
mode whirl. The sketch indicates where measurements were 
taken. The oscilloscope in orbital mode shows that shaft in 
whirl at lower rotative speeds vibrates in phase: at higher 
rotative speeds. shatt right-hand side is 180° out of phase 
from left-hand side. 
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Figure 5 
spectrum cascade of the vibrational response of a rotor indi- 
cating coexistence of seli-excited vibrations of the whip and 
second mode whirl types. 
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pear, then become “negative”, More precisely, the actual 
damping forces do not change at all, but there are emerging 
forces which act opposite to damping, first nullifying it, then 
becoming dominant. 

The “stable” rotor at rotative speeds below the threshold of 
stability rotates smoothly, and usually, due to residual unbal- 
ance, it exhibits some synchronous lateral vibrations (and 
possibly higher order harmonics). At the threshold of stahil- 
ity, the rotor starts vibrating with a subsynchronous fre- 
quency and increasing amplitude. finishing up in a limit 
cycle self-excited vibration of the whirl or whip category 
(Fig. 6h While this new subsynchronous component 
appears in the vibrational spectrum, the synchronous com- 
ponent docs not exhibit any significant change, 





Figure 6 
Oil whirl limit-cycle inception. Numbers on the orbit indicate 
consecutive rotations. Picture taken from oscilloscope display 
with exposure 1/3 second; shaft rotative speed was 3180 rpm. 


The threshold of stability can be easily predicted by a linear 
mathematical model of the rotor system (such as represented 
by hearing or seal coefficients), However, the self-excited 
vibrations can only be described by a nonlinear model. 
The nonlinear terms in the model determine the values of 
vibration amplitudes, i.e., the size of the whirl and/or whip 
limit eyele orbits. 

7. How to predict the threshold of stability 
The threshold of stability does not depend exclusively on the 
bearing/seal characteristics; it 1s definitely an entire rotor! 
bearing’ scal system property. 
For a laterally isotropic rotor supported in one rigid bearing 
and one oll-lubricated bearing (Pig. 7), the threshold of sta- 
bility, i2g7, can be approximated [Ref. 15] as: 
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where A) and A» are the partial sulfnesses of the rotor first 
lateral mode, Af is the modal mass. Ag 1s fluid film radial 
stiffness at zero eccentricity, and A is the fluid circumferen- 
lial average velocity ratio. Since the stiffMesses Ay and Aa 
act in series, their equivalent stiffness is smaller than either 
of them. Therefore, when cither Ap or As is small, their 
cquivalent stiffness becomes insignificant in comparison to 
Ki, and can be neglected. This yields the final approxima- 
tion in Equation (1). 


In the case of a rotor with a close to mid-span seal and two 
rigid supports (Fig. &), the threshold of stability is relatively 
higher than (1): 


(K, + Kom) 





Usp = i (2) 


where A), As and Af are the corresponding partial stitt- 
nesses and mass of the shaft first lateral mode. 


In both of these simple cases, the threshold of stability 19 in- 
versely proportional to the fluid circumferential average 
velocity ratio, A, and to the square reot of rotor modal mass. 
Since A decreases with shaft eccentricity inside the bearing 
or seal, the threshold of stability increases with eccentric- 
ity. This well known that when a high radial load 1s applied to 
the shaft (such as gravity for horizontal machines}, it results 
in shaft displacement to higher eccentricity inside the bear- 
ings. This stabilizes the rotor: oil whirl or whip disappear, 

Determined by the rotor actual mass, the modal mass Af in 
the denominator of the threshold of stability is also signifi- 
cant: a heavier rotor usually results in lowering the thresh- 
old of stability. This sounds contradictory to the common 
beliefs. Note, however, that we are talking about the “muss” 
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not the “weight” (as a force) and we do not specify whether 
the rotor is horizontal or vertical. It is assumed that the shaft 
rotates at the same original eccentricity, i-c., A remains the 
same for the rotor with lower or higher mass. 


Interesting also is the influence of rotor stillness on the sta- 
bility threshold. If at the rigid bearing side in the first model 
(Fie. 7) the rotor is stiff (A) higher), the threshold of stabil- 
ity increases. [f it is stiff at the fluid lubricated bearing side 
(ie., the rotor mass center is closer to the fluid bearing}, the 
stability threshold decreases, i.c.. the rotor is less stable 
(Fies. 9 and 10). 


What are self-excited vibrations? 


Vibrations in mechanical systems are usually classified into 
three categories: free vibrations, forced vibrations, and self- 
excited vibrations (Fiz. 11). The differences between them 
exist in the mechanism of energy supply to sustain the 
vibration, 
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Spectrum cascade of the rotor vibrational response during 
startup, Heavy disk clase to the rigid bearing. K] high. thresh- 
old of stability high. 
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spectrum cascade of the rotor vibrational response during 
startup. Disk close to the oil-lubricated bearing. A, low. 
threshold of stability low. 


Figure 8 
Model of a rotor!/seal system. 
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Free vibrations, forced vibrations, and self-excited vibrations 


in mechanical systems. 


Free vibrations are of the transient type and they follow an 
Instantancous energy impulse, For lightly damped systems, 
the dominant frequency of free vibrations corresponds to the 
lowest natural frequency of the system. 

In order to excite forced vibrations, the applied external 
force (thus the energy input) is usually periodic: the external 
forced vibralion frequency is equal to that of the exciting 
force. 

Note that both of these vibration types can exist in rotating or 
nonrotating structures. In general, the input force is inde- 
pendent trom rotative motion, The resultive vibrational 
mation does not affect the input force, 


selexcited vibrations are induced by a constant force, 
sustained by a constant energy supply. The system has an im- 
ternal energy transfer mechanism which delivers the encrey 
Ina periodic manner, The frequency at which the energy is 
provided usually corresponds to a natural frequency of the 
system. 


A stable system responds to an impulse force with decaying 
vibrations: an unstable system responds with increasing vi- 
bration amplitude. Theoretically for unstable linear systems, 
this amplitude grows to infinity. Practically, the amplitude 
increase is accompanied by a significant increase of system 
nonlineariues which, in turn, slows down the amplitude 
erowth. Most often, the final result is a “limit eyele” self-ex- 
cited vibration. The amplitude of the self-excited vibration 
is, therefore, limited, and determined by nonlinearities in the 
system. The frequency of self-excited vibrations is very 
close to the linear system natural frequency at the threshold 
of stability, as it is usually insensitive to nonlinearities. 
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In rotor/bearing//seal systems, the constant enerey input 1s 
represented by rotation. Fluid dynamic forces generated by 
shaft rotation in bearings and seals act as enerey convertors 
from rotation to lateral motion and cause rotor self-excited 
vibrations of the whirl’whip type. Following all classical 
features of self-excitation in mechanical systems, these self- 
excited vibrations have a frequency close to one of the sys- 
tem natural frequencies, and the amplitude is determined by 
the system nonlinearities. The most dominant and signifi- 
cant among these nonlinearities are fluid film radial stiffness 
and fluid circumferential average velocity ratio as nonlinear 
functions of shaft eccentricity. 


Example: Natural frequencies and whirl/whip 
amplitudes of the rotor supported in one rigid 
and one oil-lubricated bearing 

Ina very simplified form, the characteristic equation for nat- 
ural frequencies and stability threshold of the lincar rotor! 
bearing system illustrated in Fig. 7 is as follows (the formal 
analysis is given in | Ref. 15]): 


(Ki + Ko - Mat) [K+ Kot iD (oA) -Ki = 0 (3) 


where w is the complex eigenvalue, Ag is the linear part of 
the fluid Alm radial stiffness, D is Quid film radial damping 
and {2 1s rotative speed. This equation has three roots of 
which the real parts representing the system natural frequen- 
cies are approximately equal to the following values: 


(4) 


i) = Aft (5) 
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The natural frequency (4) corresponds to the rotor first 
lateral mode (“whip frequency 7). The natural frequency (3) 
is due to the fluid'solid interaction and is called whirl 
frequency ` 


The equation to calculate the frequency of the sell-excited 


Vibrations differs very little from equation (3): 
(AK) +o - Afar) [Kat Apt FAIA fe] 


(6) 


1 
- K =0 


where Yo = Ay =- Ko is a nonlinear part of the fluid film 
radial stiffness versus eccentricity (FG) =0, Fig. 1} and A 
is the Journal self-excited vibration amplitude (“a rotating 
eccentricity” in this case}. The self-excited vibration fre- 
quencies will differ, therefore, very little from the natural 
frequencies (4) and (3). In this example, the fluid circumter- 
ential average velocity ratio A is assumed to be constant. 
The equation resulting from (6) for journal self-excited vi- 
bration amplitude estimation is as follows: 


Mar-K 


Kr Ky “+ Fid) = Ka (7) 


A + KasA 


April 1989 


In equation (7), A and œw are journal amplitude and fre- 
quency of the self-excited vibration correspondingly. 

The whirl amplitude can be estimated when the self-excited 
frequency is w = AQ, obtained from the lingar analysis, 
Equation (7) reads, therefore, (Fig. 12): 


MR - K] 
Ag+ VLA whirt) = Ka o (8) 


Ky + KMKR? 


Note that for a rotative speed lower than the threshold of sta- 
bility, i.e., when 





AY Koka 


Mo MiKy+ Ka) 
the right-hand side function of (8) is smaller than Ag + Y, 
and the two functions do not intersect (Fig. 13). This per- 


fectly corresponds to the previously estimated stable rotor 
behavior below the threshold of stability. 









= gT 


For Q> Qer both functions of Equation (7) cross, and at the 
intersection, yield the whirl or whip amplitude (Fig. 13). 
When particular numerical values of the functions F and A 
are given, Al whi can be calculated numerically, Usually the 
nonlinearity of A at low and medium eccentricitics intro- 
duces only a slight numerical adjustment, without modity- 
ing the qualitative picture of the phenomena. 


The right-hand side function in Equation ($) increases to m- 
finity when { approaches the value VWOA)+Aa)/M4. The 


latter value corresponds to the transition between “whirl” 
and “whip” natural frequencies (5) and (4) respectively. 
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Figure 12 
Estimation of whirl or whip amplitude based on Equation (7). 
Ag = Ag + ¥(-1) is fluid redial stiffness as a function of eccen- 
tricity: cis bearing clearance. 
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Figure 13 
Estimation of the whirl and whip amplitudes (A) for various 
rotative speeds (9) and system naturel] frequencies w, based 
on Equation (7). 
From Equation (7), it follows also that when w equals the 
whip frequency (4), the right hand side function (7) tends to 
infinity. Since the bearing radial stiffness Xp tends to infin- 
ity when the shaft approaches the bearing or seal wall, the 
journal whip amplitude also approaches the bearing wall. 
Thus the journal amplitude is limited by the bearing/seal 
clearance, c, and A pip = C. 
More important in the whip case is, however, to analyze the 
value of the rotor amplitude at mid-span, as the rotor 
vibrates at its resonant frequency. For both whirl and whip, 
this amplitude can be calculated from the following 
equation: 
Ka 
Arao = A journal DaannO a (9) 
VIK AK -MW + Dear 
with the relative phase between journal and rotor vibrations 
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Š = arctan = al 
Where Dg is rotor modal damping, 


Por whirl, ic. when œw = AG: 
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The phase @ is close to zero, The whirl mode has a “coni- 
cal” shape (Fiz. 14). 


For whip, when a = WA) +K) OM and A journgs = E, M 


follows from (9) and {10} that 





Ka 

Aroor = cek + (12) 
Dy 

g= 90°, (13) 
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thus, the rotor amplitude can be very high as it 1s controlled 
only by the rotor damping Dy. Note that the bearing damping 
D has oo influence on the whip amplitude. The rotor vibrat- 
ing at whip lrequeney is in resonance and, therefore, is in a 
highly dangerous condition. This condition is similar to pro- 
longed rotor operation at the first balance resonance speed. 
But this time this jeopardizing condition is rotative speed 
independent. “Passing through,” like the first balance reso- 
nance critical speed, is not possible. There exist. however, 
some other measures to cormeet the problem. 
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Higher eccentricities also lower the A ratio which further 
improves stability, It requires an increase in the radial force 
applied to the shaft. It can be created cither by a “friendly” 
misalignment, and/or by a heavier shaft in horizontal 
machines and/or by main flow distortion in fluid handling 
machines, so that the flow generates a significant radial 
force component, 


An increase of Muid film radial stiffness, Aig, can also be 
achieved by increasing fluid pressure. Hydrostatic bearings 
(or more properly called “externally pressurized” bearings} 
are well recognized for their high stability features, 

These cures which prevent rotor whirl and whip are 
lustrated in Figs. 15 to 18. The relationships are deduced 
trom Equation (7). 
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Whirl and whip modes of the rotor supported in one rigid and 
one oil-lubricated bearing. Figure 15 

Instability cure: threshold of stability increases for shafts 
with fluid anti-swirl injections: A gra is the fluid average veloc- 
ity ratio due to anti-swirl. 


10. Whirl and whip cures 


Analyzing the threshold of stability (1) and (2) and whorl! 
whip amplitude evaluation (Equation (7) and Fig. 13), it ts 


casily scen that there exist two major possibilities to im- a pad = tfa 

prove the rotor stability, independently of rotor parameters: aa 7 | 

Decrease of A and increase of Ay. m Hi eee 

A decrease of the fluid circumferential average velocity ra- a a 5 fa | | 

to, A, can be achieved by a perturbation of the circumteren- = Mt H> 

tal flow regular pattern, or prevention of the cireumlerential = a 7 l | T 

flow occurrence in the bearing for seal. In bearings, this ae E 

can be successfully accomplished by noncircular geometries BESSLNG sem 

ceg., elliptical, two or three lobe) andor the presence of CENTERLINE = = ‘ORIGINAL THRESHOLD 


moving parts inside bearings (iltine pads, floating rings). In HE a ha 






such bearings, the circumferential flow ls significantly x / is | 

Teduccel in comparison to circular bearings. i* F | a HIGHER ECCENTRICITY 
. SETT = fal OPERATION 

In seals, a decrease of A can be accomplished by injections Agy M 


ies ee aie ) | Aa Ao 
of the Huid tangentially in the direction opposite to rotation 

In order to perturb the rotation-generated circumferential a ee 
flow. This technique is known as “anti-swirl” control [Ref. 0 es | z 
[8,19] and is currently widely used in compressors, turbines, SHAFT ECCENTRICITY 
and pumps. Figure 16 
Instability cure: threshold of stability increases for shafts ro- 
tating at higher eccentricity: Aq is the fluid average velocity 
ratio at higher eccentricity. 


Moving the shalt to a higher eccentricity affects the bearing’ 
seal racial stiffness Ay. At higher eceentriciuly, the radial 
stiffness is much larger and this improves rotor stability, 
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Figure 17 
Improvements of rotor stability with externally pressurized 
bearings; Ap reaches a higher value than the threshold of ster- 
bility. 
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Figure 18 
Instability cure: whirl'‘whip occurs when the shaft is concen- 
tric in the bearing or seal. Improvement of rotor stability is 
achieved when shaft operates at higher eccentricity: Func- 
tion Ko(Mow*-K)) Í (K,+Ke-Mw*) from equation (7) appears 
below the new Ko and does not cross with Xa: Ne whirl! whip 
exists. 


ll. Closing remarks 
This article presents new results regarding flow -induced 
instabilities of roter/bearing/seal systems. Experimentally 
identified models of fluid forces allow tor better understand- 
ing of vibrational phenomena occurring in rotating machines 
and lead to rational prevention of undesirable (low-induced 
self-excited rotor vibrations. 
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Notation 
Amplitude of rotor self-excited vibration of whirl 
or whip type. 
Amplitude of journal whirl or whip, 
Amplitude of rotor mid-span whirl or whip. 
Amplitude of journal whip. 
Amplitude of journal whirl. 
Bearing or seal radial clearance. 
Bearing or seal radial damping coefficient. 
Rotor first lateral mode damping coefficient. 
Y-i 
Bearing or seal fluid film radial stiffness, as func- 
tion of shaft cecentricily. 
Bearing or seal fluid film radial stiffness at zero 
eccentricity, 
Rotor first lateral mode partial stiffnesses. 
Rotor first lateral mode generalized mass, 
Time. 
Rotor lateral detlection (x — horizontal, y — 
vertical, 


Relative phase between rotor and journal self-ex- 
cited vibrations. 


Fluid circumferential average velocity ratio. 


Fluid circumferential average velocity ratio for 
concentric shaft, 


Anti-swirl-related fluid circumferential average 
velocity ratio. 

Bearing or seal Awd film radial stiffness, a nonlin- 
ear function of shalt eccentricity. 


Rotor’ bearing /seal system natural frequency, also 
frequency of self-excited vibrations. 


Rolative speed. 
Threshold of stability, 
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